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Abstract. The structures observed in the sub-threshold neutron- induced fission of 
^^^Th were investigated employing a recent developed model. Theoretical single- 
particle excitations of a phenomenological two-humped barrier are determined by 
solving a system of coupled differential equations for the motion along the optimal 
fission path. A rather good agreement with experimental data was obtained using a 
small number of independent parameters. It is predicted that the structure at 1.4 and 
1.6 MeV is mainly dominated by spin 3/2 partial cross-section with small admixture 
of spin 1/2, while the structure at 1.7 MeV is given by a large partial cross section of 
spin 5/2. 
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1. Introduction 

The measured neutron-induced fission cross-section behavior of nuclei in the thorium 
region represented a challenge for nuclear physicists concerning the shape of the potential 
energy surface. The experimental data suggested the existence of a triple humped 
barrier. The neutron-induced cross- sect ions of 230,232r|n]^ exhibit multiple fine structures 
[H [21 [3] superimposed on a gross structure of the threshold cross-section. If the fine 
structure is interpreted as a serie of rotational states constructed on a /3-vibrational 
state produced in some well of the deformation energy, it is straightforward to postulate 
the existence of a triple-humped barrier. The spacing between the members of the band 
is so small that it is consistent only with a prolate deformation that reaches the vicinity 
of the second-barrier top. The analysis of Ref. [1] indicates that an intermediate state 
nucleus must exist at a deformation considerably larger that the normal value. A ternary 
minimum obtained theoretically in the potential energy surface of ^^'^Po [4] made this 
hypothesis credible. Therefore, a shallow minimum was assumed at this deformation 
to create a new /^-vibrational state. Angular distribution analysis [H [6] confirmed the 
existence of the triple well. Up to now, the assumption of a triple-humped barrier 
seems to be the best interpretation for the fine structure of intermediate cross-section 
resonances |7]. 

On the other hand, our analysis explores a different way to consider the cross-section 
resonant structure phenomenon by quantifying the dynamical single-particle effects 
associated to vibrational resonances produced in the second well [8]. Our exploratory 
investigation showed that the ^^°Th neutron-induced fission threshold resonant structure 
can be explained [9] by rearrangements of single-particle orbitals on the way from the 
initial configuration of the compound nucleus up to scission. This resonant structure 
depends also on the dynamics of the process. 

Sect. [2] will provide a general description of the formalism intended for the 
evaluation of single-particle excitations, while results concerning the intermediate 
structure of the fission cross-section will be extensively presented in Sect. [H Comments 
are made in Sect. HI 

2. Single-particle excitations 

In most usual theoretical treatments of nuclear fission, the whole nuclear system is 
characterized by some collective coordinates associated with some degrees of freedom 
that determine approximately the behavior of many other intrinsic variables. The 
basic ingredient in such an analysis is a shape parametrization that depends on several 
macroscopic degrees of freedom. The generalized coordinates associated to these degrees 
of freedom vary in time leading to a split of the nuclear system in two separated 
fragments. A microscopic potential must be constructed, to be consistent with this 
nuclear shape parametrization. The three important degrees of freedom encountered in 
fission, that is, elongation, necking and mass- asymmetry, must be taken into account. 
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By solving the Schrodinger equation for a reasonable mean field potential associated to 
the nuclear shape parametrization, the single-particle energies are determined. In the 
case of odd-nucleon systems, the potential barrier must be increased with an excitation 
associated to the unpaired nucleon. The amount of which the barrier is increased can 
be estimated within the specialization energy p^. This quantity can be interpreted as 
the excess of the energy of the unpaired nucleon with a given spin over the energy of 
the same spin nucleon state of lowest energy. 

In the present work, an axial-symmetric nuclear parametrization is obtained by 
smoothly joining two intersected spheres of different radii Ri and R2 with a neck surface 
generated by the rotation of a circle of radius Rs around the symmetry clXlS, clS displayed 
in Fig. [H The surface equation is given in cylindrical coordinates: 



where Zd and Zc2 define the region of the necking. The meaning of the geometrical 
symbols that depends on the shape parametrization can be understood inspecting Fig. 
[H This parametrization allows to characterize a single nucleus or two separated nuclei. 
Throughout the paper, the subscripts 0, 1, and 2 indicate the parent, the heavy and 
light fragments, respectively. If 5"=!, the shapes are necked in the median surface 
characterizing scission shapes and if 5*=-! the shapes are swollen characterizing the 
ground-state and saddle points. The macroscopic parameters used in the following 
are denoted R = Z2 — zi (elongation), C = S/R3 (necking) and 77 = R1/R2 (mass- 
asymmetry). For large distances between the two nascent fragments, the configuration 
given by two separated spheres is reached. 

For instance, to minimize the action integral [H] it is very difficult to treat the three 
independent generalized coordinates in the same time. Some simplifying assumptions 
must be introduced. As mentioned also in Ref. [12], microscopic approaches to fission 
[131 E] established that the second saddle point is asymmetrical with a value compatible 
with the observed mass ratio. In the same time, in the region of the second barrier, the 
mass-asymmetry component of the inertia tensor is very large [I5]. So, the variations 
of the mass-asymmetry coordinate are hindered in this region. On another hand, 
for elongation smaller than that of the second well deformation, the mass- asymmetry 
component of the inertia is much lower. Therefore, the mass-asymmetry coordinate can 
be modified without enhancing too much the value of the action integral. Moreover, the 
deformation energy is less sensitive to variations of the mass-asymmetry coordinate in 
the region of compact shapes. As in Ref. |15|, this observation allows us to reduce the 
number of parameters in order to rend our problem tractable. Therefore, the evolution 
of the mass asymmetry generalized coordinate will be a priori fixed in the following. It 
is assumed that the ratio R1/R2 varies linearly from unity (first barrier top) to the value 
associated with the final mass partition (second barrier top). The mass asymmetry in 
the outer barrier region is deduced by considering that the volume occupied by the light 




(1) 
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Figure 1. Nuclear shape parametrization. zi, and Z3 are the positions of the 
centers of circles of radii i?2 characterizing the two nascent fragments, and of i?3 
determining the neck, respectively. If s = 1, the shape is necked, otherwise the shape 
is swollen in the median surface. The distance between the two centers z\ and zi 
determines the elongation R. 



fragment equals the final one. 

The deformation energy of the nuclear system is the sum between the liquid drop 
energy and the shell effects, including pairing corrections. The macroscopic energy is 
obtained in the framework of the Yukawa-plus-exponential model extended for binary 
systems with different charge densities The Strutinsky prescriptions [17] were 

computed on the basis of the Superasymmetric Two Center Shell Model (STCSM) 
p!8t [T9]. For one of the most probable partition ^^^Th — > ^'^Sr -|-^^^Te, the deformation 
energy as function of C and R is plotted in Fig. [2l 

The theoretical study of binary disintegration processes is limited by the difficulties 
encountered in the calculation of single-particle levels for very deformed one-center 
potentials. On one hand, central oscillator potentials are not able to describe in a 
correct manner the shapes for the passage of one nucleus to two separated nuclei without 
including a large number of multipole deformation parameters and, on the other hand, 
for very large prolate deformations the sum of single-particle energies reaches an infinite 
value, as evidenced within the deformed oscillator model. These difficulties are surpassed 
by considering that the mean field is generated by nucleons moving in a double center 
potential. This kind of models allows to describe scission configurations within a small 
number of degrees of freedom. A more realistic version of the two-center shell model was 




zcl 
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realized recently [IS] and it is used to generate the single-particle energy evolutions from 
the ground-state up to the formation of two separated fragments. The nuclear shape 
parametrization being characterized by an axial symmetry, the good quantum numbers 
are the projection of the spin Q. As in the Nilsson model, the single-particle energies 
depend on two interaction constants k and rj, related to the spin-orbit operator and to 
the squared orbital momentum correction, respectively. These constants are determined 
in order to reproduce the ground-state properties pO] . 

In order to determine the single-particle excitations, it is not sufficient to have a 
model for the intrinsic nuclear levels, but is necessary to perform a full calculation of the 
trajectory of the decaying system in the configuration space. The shape of the fission 
barrier can be obtained if the trajectory of the nuclear system in our three-dimensional 
configuration space is obtained, starting with the ground-state of the compound nucleus 
and reaching the exit from the second barrier. This trajectory can be obtained by 
minimizing numerically the action functional that gives the quantum penetrability: 



in the semi-classical Wentzel-Kramers-Brillouin approximation [T7]. The two turning 
points Ri and Rf denote the elongations that characterize the ground-state and the 
exit point of the barrier, respectively. Here V{R, C, rf) is the deformation energy and 
M{R,C,r], 1^, 1^) is the effective mass along the trajectory. The inertia is computed 
in the frame of the Werner- Wheeler approximation [21j, that means, the flow of the 
fluid is idealized as non-rotational, non-viscous and hydrodynamic. Using the minimal 
action principle, in general, the nuclear system does not follow a path characterized by 
minimal values of the deformation energy, so that the trajectory does not interpolate 
barrier saddle points values. 

Having in mind the assumption imposed for the variation of the mass- asymmetry, 
the action integral must be minimized in a two-dimensional space spanned by C and 
R. The flrst turning point Ri is flxed but the second Rf lies on the equipotential line 
that characterizes the exit from the outer barrier, that is i?j is a function of C. A 
simple numerical method is used to flnd the paths characterized by different values 
of Rf, associated with local minimums. For that purpose, the function C = f{R) is 
approximated with a spline function of n variables Cj {j = 1, n) in flxed mesh points 
Rj located in the interval [Ri,Rf] along the elongation axis. A numerical expression 
for the WKB functional ([2]) that depends only on the parameters Cj is obtained. This 
expression is minimized numerically. For every value of Rf a. local minimum is obtained. 
The best values are retained. The trajectory is displayed on Fig. O This dynamical 
trajectory starts from the ground-state, reaches the region of the second well and the 
slope changes suddenly to penetrate the outer barrier. Between the flrst and second 
well, the macroscopic coordinate C is less than 0, that is the shapes are swollen in the 
median region. Penetrating the second well, the shapes become necked. The theoretical 
potential barrier obtained along the minimal action path is plotted in Fig. [31 The 
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Figure 2. Deformation energy in MeV for the partition ^^■^Tli^^^Sr+^'^^Te. C 
represents the curvature of the neck and R the distance between the centers of the 
fragments. Positive values of C characterize necked-in shapes. The mass-asymmetry 
is varied linearly with R from a value ry(i? ~ 5 fm)=0 (close to the ground-state of 
the compound nucleus) to the final value rj — A1/A2 (in the vicinity of the top of the 
second barrier). The step between two equipotential lines is 1 MeV. Several values of 
the deformation energy are marked on the plot. The dynamic trajectory is represented 
with a thick line that starts in the first well, penetrates the first barrier, attains the 
second well and finally tunnels the second barrier towards scission. 

height of the outer barrier is very large, therefore some corrections are required in order 
to obtain reahstic values of the fission cross-section. This is the main reason that leads 
to use a phenomenological barrier in calculating the cross section. The first well is 
located at approximately R=5.5 fm and identifies the fundamental state. In Fig. HI the 
nuclear shapes of the extreme values of the barrier are displayed. 

Using the STCSM the neutron diagram is computed along the minimal action 
trajectory, as displayed in Fig. [51 Up to i? ~ 5.5 fm the nuclear system is considered 
reflection symmetric. From the ground-state (located at approximately -R=5.5 MeV) 
up to scission, the system loses the reflection symmetry to reach the final partition 
233rp]^_j^98gj,_j_i35'jg_ thesc circumstanccs, the parity is no longer a good quantum 
number, the levels being characterized only by the spin projection Q as good quantum 
numbers. The Nilsson coefficients of the orbital momentum operators (k=0.063 and 
77=0.8) were determined to reproduce as better as possible the experimental sequence 
of the first excited levels in ^^^Th. The first single-particle excited states are retrieved: 
an 1^ state (fundamental level) emerging from 2g9/2 followed by a |^ one. 

To determine the cross section, several single-particle levels are selected that lie as 
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Figure 3. Theoretical dynamical barrier calculated along the minimal action 
trajectory as function of the elongation R. 



(a) (b) (c) 




(d) (e) 

Figure 4. The shapes obtained along the minimal action trajectory, (a) The ground 
state with elongation i?=5.8 fm and necking coordinate C=-0.053 fm~^. (b) The region 
of the first barrier with i?=10.57 fm and C=-0.04 fm~^. (c) The region of the second 
well with i?=13.69 fm and C=-0.0508 fm~^. (d) The region of the second barrier with 
i?=15.139 fm and C=-0.008 fm^^. (e) The region of the exit from the barrier with 
R=17 fm and C=0.085 fm-^. 



close as possible to the Fermi energy region. These levels give the major contribution in 
the strength of the fission channel due to their low excitation energy and the large 
amount of macroscopic kinetic energy available for fission. Concerning the Q=l/2 
workspace, 8 selected levels, Ei up Eg are extracted separately in the left panel of 
Fig. [6] as an example. The last occupied level is denoted Ep. The diabatic levels 
of the subspace Q = 3/2 are displayed in the right panel of the same figure. In the 
following, for simplicity, the discussion will be restricted only for the subspace Q = ^. 




Figure 5. Neutron level scheme as function of the elongation. At elongation zero, 
the shape parametrization gives a spherical nucleus and the spectroscopic notations 
are available. For low values of the deformations, the system behaves as a Nilsson 
level scheme. Asymptotically (R — > oo) the two diagram of the two formed fragments 
are superimposed. In the adiabatic representation, the last occupied level is displayed 
with a thick dashed line. The 8 selected levels with f2=l/2 are represented with full 
thick line, the 5 levels with $1=3/2 arc plotted with dot-dashed thick lines, the 4 times 
il=5/2 and 3 times 7/2 levels are marked with dotted lines (smaller distance between 
points for Cl=5/2). The ground-state of the compound nucleus is indicated with an 
arrow. 

For Q = |, |, |, the same procedure as in the case of Q = | will be used. 

A first behavior can be noticed. The nucleon located on the adiabatic level emerging 
from El reaches a very unfavorable energy configuration after the scission. In the 
fundamental state, this unpaired nucleon is located on the fundamental level but arrives, 
in the adiabatic representation, at several MeV under the last occupied level (the €3 
diabatic level). So, if the nucleon is initially on the ground-state, it must follow a 
diabatic energy path to arrive in a most favorable energy configuration, that is close to 
the last occupied level (in one of the diabatic states Cq, 65, ey or £2). So, adiabatically, 
the fission strength for states with spin 1/2 is not favored. This effect is a direct 
consequence of the rearrangement of low spin orbitals during the disintegration. The 
number of levels with = | in the two nascent fragments that are under the energy of 
the last occupied level is always larger that the same number in the compound nucleus. 
So, fl=l/2 orbitals with larger energies of the parent must decrease in energy to fill 
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Figure 6. (a) The 8 selected levels with fl—1/2. The last occupied level in the 
adiabatic representation is denoted Ep and is represented with a dashed line. The 
avoided level crossing regions are numbered and the diabatic levels identified. In the 
ground-state configuration, the ei level (emerging from Ei) is superimposed to Ep- 
At i? ~ 20 fm, Ep is located between eg and £5, while the adiabatic level emerging 
from El dropped to £3. (b) As in plot (a) for the 5 levels with fi=3/2. With thin dot 
dashed lines the 4 f7=5/2 adiabatic levels are also displayed. 



the levels located under the Fermi energy of the two fragments. This aspect somewhat 
hinders the possibility to fission through ^ = \ channels. The next step is to study the 
energy paths followed by the unpaired nucleon in the diagram. 

The realistic two-center diagram presented before provides an instrument to study 
the role of individual orbitals during the disintegration process in a similar way as 
the study of nucleus-nucleus collisions [221 [23] or the alpha- and cluster-decays [2H [25] . 
Levels with same quantum numbers associated to some symmetries of the system cannot 
cross during the disintegration process and exhibit avoided level crossing. In our case, 
due to the axial-symmetry of the system, the good quantum numbers are the projection 
of the spin Q. The point of nearest approach between two levels of same fl define 
an avoided level crossing region. If the internuclear distance varies, the transition 
probability of a nucleon between two adiabatic levels is strongly enhanced in the avoided 
level crossing region. This promotion mechanism is generically called the Landau-Zener 
effect. 

Concerning the 8 single-particle adiabatic levels {Ei, Eg) belonging to the Q=l/2 
workspace, the first step is to find the avoided level crossing regions. The avoided 
crossing regions can be obtained by plotting the energy differences between these 
adiabatic levels as in Fig. [71 Each pertinent avoided crossing is identified and numbered. 
The avoided crossings that have a chance to be located along the diabatic single-particle 
energy paths emerging from the lower levels Ei, E2 and E^ are considered pertinent. 
Due to their low initial excitation energy, the transitions through these levels carry the 
major part of the fission strength. That property allows to restrict our calculations only 
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Figure 7. Differences between tlie selected adiabatic levels. The avoided level crossing 
regions that appears between the adiabatic energies emerging from the initial states 
El, ...E4 are numbered as in Fig. [6l 



for an initial condition in which the occupation probabihty of one of these levels is one. 
The next step is to determine the probability of realization of each diabatic energy path 
emerging from these levels. Concerning the f2=3/2 subspace, the analysis is realized in 
the same manner, for initial conditions restricted to the first 3 low energy levels. 

Assuming an n-state approximation, the wave function of the unpaired nucleon can 
be formally expanded [26] in a basis of n diabatic wave functions 0i(r, R) as 



^'(r, R,t) = J2 Ciit)(l)i{r, R) exp (^-^ ^ e^dt 



(3) 



where the matrix elements with the diabatic states are abbreviated as follows 

eij =<(f)i\H\ (f)j > (4) 

where H is the STCSM Hamiltonian and q are amplitudes. Inserting \E' in the time- 
dependent Schrodinger equation. 



H 







(5) 
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Figure 8. The occupation probabilities of the diabatic levels ei, ...eg as function of 
the distance between the centers of the fragments. The same numbers as in Fig. [7| 
are used to identify the avoided level crossing regions. The internuclear velocity is 
3.5x10'' m/s, that leads to a reasonable reaction time (time to penetrate the barrier) 
of approximately 5x 10~^^ s. This example is constructed for an initial condition 
while p^. — {i ^ 1). The occupation probabilities vary in the avoided level crossing 
regions. 

the following system of coupled equations is obtained: 

^^ = -f^Yl Cj-eij exp f - - / {ejj - eii)dt j (6) 

To solve this system, the internuclear velocity the diabatic energies and the 
interaction matrix elements must be known. Apart the relative velocity, the other 
ingredients are supplied by the STCSM. The diabatic states are constructed by using 
spline interpolations in the level crossing regions. The interaction matrix elements 
eij between the diabatic states is a measure of the difference between adiabatic and 
diabatic energies. The occupation probability of each adiabatic level as function of R is 
now obtained p^. =| c, p. For the unpaired neutron initially located in the fundamental 
state El, the system ([6]) is solved within the boundary condition ci=l and Ci=0 for i ^ 1. 
The occupation probabilities of each diabatic level plotted in Fig. [6] are represented in 
Fig. [H Within the selected levels and avoided level crossings, 40 different energy paths of 
the unpaired neutron can be obtained as indicated in Table [H Here, an approximation 
is made by considering that the points of the avoided level crossings 1 and 2 form a 
single avoided level region. Otherwise, the number of paths gets doubled. Each path 
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Figure 9. The fl—1/2 phenomenological Vph barriers with excitations obtained in the 
frame of the hybrid model emerging from the fundamental level Ei . 

represents an excitation of the nuclear system. The probabihty of each path can be 
estimated. For example, it can be deduced from Fig. [8] that the path ei — 2 — ei — 3 — €4 
carries about 0.5 of the probability. The line between letters and digits connects diabatic 
levels and avoided level crossing regions. A strong mixing is produced in the region 4, 
that leads us to conclude that the paths ei — 2 — ei — 3 — £4 — 4 — £3 (no. 1 in table [I]) 
and ci — 2 — ci— 3 — £4 — 4 — £4 — 7 carry each of them about 0.25 probability. Finally, 
it can be considered that the path ei — 2 — ei — 3 — 64 — 4 — £4 — 7 — ei (no. 2) has 
about 0.05 probability of realization while ei — 2 — ei — 3 — £4 — 4 — e4 — 7 — e4 (no. 3) 
remains with 0.2. The other probabilities are estimated in the same manner. The same 
procedure is repeated for the case when the unpaired nucleon is initially located on the 
other selected levels. 

The excitations of the barriers due to one diabatic path k is given by the 
specialization energy. Considering that the fundamental barrier corresponds to the 
nucleon at the Fermi energy, the excitation i?^ as function of R is 

E.kiR) = V(efc(i?)-A(/?))2 + A2(i?) - A{R) (7) 

in the frame of the superfluid model. Here, is the single-particle energy of the path k, 
X is the Fermi energy and A the gap. These excitations are added to the fundamental 
barrier. These quantities have the same meaning as the so-called transition bandheads 
found in the literature. 
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3. Cross-section 

The partial fission cross section ctj for a spin J of tlie compound nucleus and excitation 
energy E* is obtained within a statistical principle: 

where the ratio in the right-hand side is the probability that the system decays through 
fission. It is given by a ratio between energy widths for fission (subscript /), neutron 
emission (subscript n) and 7-de-excitation (subscript 7). The subscript T addresses 
the total transmission in the fission channel including absorption in the second well. 
The neutron transmission was computed for a squared complex potential [27] in order 
to evaluate the compound nucleus cross section. To determine the participation of 
different Q excitations in the fission channel for a given spin J of the compound nucleus, 
a unfolding procedure in term of Clebsh-Gordon coefficients is used 
r ( T T^*') — 1 \^Lm <jLno\jn>'^ 

X /o^*"^^ 2Tf{E, L, n)p{n, E* -E- EL)dE ^ ' 

where a normalization coefficient is used: 

c = J2J2< JL^^ I >^ (10) 

L=o n 

and the condition J = L+Q is imposed. Here p represents the density of states, El is the 
the rotation energy in the fundamental state of the compound nucleus with an angular 
momentum L and Lm is the maximum orbital momentum taken into consideration. 
This formula can be obtained easily by simplifying the model underlined in Ref. [8], 
that is, neglecting the additional collective excitations as gamma, sloshing or bending 
vibrations. Analog formulas can be obtained for the 7 and neutron energy widths as 
detailed in Ref [8]. 

In the fission channel, the spin Q density of states can be shared as function of the 
excitation energy between a discrete component and a continuum one: 

""'"'-l Pac(n,E), E>E, 

where pcc is the statistical Gilbert and Cameron approximation and e^^j {i = l,n) are 
the set of diabatic single particle energies that are taken into consideration for a spin 
projection Vt. So that, the transmission in the fission channel can be decomposed as 
follows: 

^E*-Er. j.^^^^ n)p{n, E* ~E- EL)dE 

= E^TfiE*-EL-en,) (12) 
+ Tf{E, L, Q)p{Q, E*-E- EL)dE 

The sum over i takes into account all the transmissions for diabatic levels with spin 
J = VL + L located in the energy interval [0,i?o]- In this context, the transmission 
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Tf{E* — El — en,i) means a weighted sum of the transmissions of all available diabatic 
energy paths emerging from the level e^^j. 

The microscopic model used to compute the theoretical barrier is subject to some 
limitations as described in Ref. [19]. It is not possible to obtain pertinent values of the 
heights of the barriers. In this circumstances, it is necessary to use a phenomenological 
barrier. A phenomenological barrier is conventionally simulated as a function of a 
dimensionless parameter (3, that characterizes a deformation, within three smoothed 
joined parabolas [28]. In our work, an imaginary component of the potential is added 
between the turning points of the second well, in order to simulate the damping 
due to gamma and neutron emission. The additional excitations are considered as 
specialization energies and are added to the phenomenological barrier. This operation 
is achieved in the simplest possible way, by realizing a linear interpolation based on a 
correspondence between the elongation R and the dimensionless parameter (3 in some 
points. The correspondence was chosen for the two minimums, the two heights and the 
exit point. The hybrid model emerges. New barriers are constructed as displayed in 
Fig. [9l When only the collective rotations are taken into account, the heights of the 
barriers and that of the second well are modified with a quantity 

AE,^^^^±^^-E, (13) 

where Ij is the moment of inertia, j labels one of the two heights or the second well. 
The decoupling parameter is neglected. The moment of inertia is computed simply with 
the formula Ij = fiR^ where /z is the reduced mass and Rj is the theoretical elongation 
obtained at the extreme point j. The quantity 

addresses the fundamental state of the compound nucleus. The previous formulas 
represents an improvement of the formalism found in Ref. [S]. 

A large number of excited states are obtained that are characterized by the 
projection Q and the angular momentum L. The transmission is calculated numerically 
by approximating the shape of the excited barrier within 500 constant potential steps 
using the numerical recipe found in Ref. [22] • A search of the heights and of the widths 
of the phenomenological barrier is realized in order to reproduce as well as possible the 
experimental fission cross-section threshold structure. A behavior that agree satisfactory 
with the experimental data is obtained. The heights of the inner phenomenological 
barrier, the second well and the outer barrier are 6.81, 4.83 and 6.61 MeV, respectively. 
In the same order, the widths are 1.2, 0.4 and 1.1 MeV. The theoretical cross-section is 
represented in Fig. [10] and compared with experimental data and evaluations. 

The evaluations succeeded to reproduce better the experimental data. In general 
many parameters are taken into account to evaluate a cross-section in terms of Bohr- 
channels. For example, in evaluations phenomenological level densities functions 
appropriately matched to the available experimental structure data at low excitation 



Threshold Resonant Structure of the Th Neutron- Induced Fission Cross Section 15 




Figure 10. Thick full-line, neutron-induced cross section for ^^^Tli as function of 
the neutron incident energy En calculated within the hybrid model. Points are 
experimental data. Thick dashed line represents the ENDF/B-IV evaluation while 
the thick dot-dashed line is the JENDL-3.3 one [3D]. Experimental data are from Ref. 
[51[31[3T]. A thin full line gives the partial cross section of spin 1/2, a dashed line is 
for the spin 3/2, the dot-dashed one for 5/2 and the dotted line for 7/2. 



energies are used. Multiplication factors are also applied to level density functions to 
account for enhancements in the fission transition state densities at each fission barrier. 
It is a common practice to describe the cross-section as the sum of excitations for 
discrete levels constructed to fit the resonance. In other words, the evaluation takes 
into account many other parameters to fit the experimental data apart the heights and 
the widths of the phenomenological barrier. In the work presented in this paper, no 
adjustments are made to improve the agreement, the simulations being based only on 
the phenomenological barrier parameters and the internuclear velocity. 

Our simulations evidence an oscillatory behavior of the cross-section close to 1.4 
MeV. This aspect is in agreement with the experimental data given in Ref. |32j. The 
experimental data combined with theoretical arguments estimate a ratio 2:1 between 
the partial cross section of spin 3/2 and 1/2, respectively. The model shows that the 
partial cross section for the spin 3/2 is responsible for the oscillations of the cross section 
at these energies. Experimentally, the peak at 1.6 MeV is explained entirely by a partial 
cross section of spin 3/2 with a small 5/2 component. In our plot a strong 3/2 component 
is present with small admixture of 1/2 and 5/2 partial cross sections. A discrepancy 
is obtained for the 1.7 MeV structure. The experiment evidences the existence of a 
mixing between 3/2 and 5/2 components while our model predicts a large 5/2 partial 
cross section followed by the 3/2 and 1/2 components. 

In Fig. [11] the cross-section is plotted on an extended scale. It can be observed 
that the theoretical results exhibits an oscillatory behavior in the low energy region, up 
to 1.2 MeV, around the smooth variation of the experimental data. In the panel (b), 
the transmissions computed for the barriers with different calculated excitations are 
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E„ (MeV) 

Figure 11. (a) Same as Fig. [10] in an extended logarithmic scale along the y-axis. 
(b) L—0 fission transmissions for different barriers as function of the neutron energy 
En- The transmissions for il=l/2 excitations are plotted with full lines, those for 3/2 
with a dashed line, those for 5/2 with dot-dashed and 5/2 with a dotted line. 

displayed. The oscillatory behavior is due to a large number of resonances associated 
to the different excited barriers. 

4. Summary and discussion 

The scope of the present work is to understand the mechanism for the formation of 
the fission cross section structure and of the high number of resonances by appealing 
essentially to dynamical single-particle effects associated to /5-vibration in the second 
well. The number of free parameters is kept as minimal as possible (six parameters 
that characterize the phenomenological barrier and one parameter for the internuclear 
velocity) to show evidence of the physics of the problem. 

Theoretical excitations and their associated probabilities were determined for a 
given partition in the isotopic distribution of fragments. These excitations are added to 
a phenomenological barrier in the framework of the hybrid model. After a suitable search 
of the parameters of the double humped phenomenological barrier, the cross section is 
computed. The results give a rather good qualitative agreement with experimental data. 
It is evidenced that the structure at 1.4 and 1.6 MeV is mainly dominated by spin 3/2 
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partial cross-section with small admixture of spin 1/2, while the structure at 1.7 MeV 
is given by a large partial cross section of spin 5/2. 

In this exploratory analysis, only one partition for the fission fragments is taken 
into consideration. For other partitions in the same mass region, it is expected that the 
level scheme changes slightly leading to a small shift in the energy of the resonances. By 
taking into account several partitions in the same mass-region and folding their yields 
it it possible to obtain broader resonances as experimentally observed. 

In this context, it will be interesting to explore experimentally if the isotopic 
fragment distribution in the fission process changes in the energetic region covered by a 
resonance, showing a preference for several mass partitions. If such a phenomenon can 
be experimental evidenced, that will represent a strong experimental support for our 
model because the statistical theories don't include ingredients related to this aspect. 

The model can be further improved. Up to now, only the radial coupling was used 
to explain the intermediate structure of the cross-section. It is possible to have better 
results by taking into account the Coriolis mixing and the residual interactions by using 
evolved forms for the system of coupled equations that describes the microscopic motion 
[251 [33]. 

Other models succeed to reproduce better the experimental data [31] using an 
extensive number of free parameters: 10 variables for the heights and widths of the triple 
humped phenomenological barrier plus 5 times 16 variables for the transition bandheads 
constructed on different intrinsic excitations (with a significance of excitations given by 
single-particle energies). Despite the overall excellent agreement on a very large neutron 
energy region, this treatment, generally used in evaluations, takes into account a peculiar 
behavior for the single-particle excitation energies. The levels that characterize the 
transition bandheads never intersect. The first |^ level have practically the same value 
(having as reference the fundamental state) during the penetration of the barrier. This 
behavior, as remarked previously, cannot be expected. Moreover, the statistical models 
consider that the population of each fundamental transition band is essentially one. 
The formalism presented in the Sect. [2] indicates that such a behavior is physically not 
reasonable. 

The present investigation shows that the resonant structure of the fission cross 

section can be explained by the existence of many barriers associated to single-particle 

excitations. So, it is possible that the complex structure in the fission cross section 

is due to rearrangement of orbitals and the djTiamic of the process, beginning from 

the initial state of the compound nucleus and terminating at the scission. A large 

number of different excited barriers are formed leading to a large number of vibrational 

resonances in the second well. These resonances carry information about the structure 

of the nucleus at hyperdeformations and the dynamics. The model presented in this 

work represents an alternative to the actual statistical models and may determine a 

competitive way to consider the fission process. 
[1] G.D. James, J.E. Lynn and L.G. Earwaker, Nucl. Phys. A 189, 225 (1972). 
[2] J. Blons, C. Mazur, D. Paya, M. Ribrag and H. Weigmann, Phys. Rev. Lett. 41, 1282 (1978). 
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Table 1. Energy paths open for the first Ei 0=1/2 level 
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